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After more than three decades the fractional quantum Hall effect still poses challenges to con¬ 
temporary physics. Recent experiments point toward a fractal scenario for the Hall resistivity as a 
function of the magnetic field. Here, we consider the so-called thin-torus limit of the Hamiltonian 
describing interacting electrons in a strong magnetic field, restricted to the lowest Landau level, and 
we show that it can be mapped onto a one-dimensional lattice gas with repulsive interactions, with 
the magnetic field playing the role of a chemical potential. The statistical mechanics of such models 
leads to interpret the sequence of Hall plateaux as a fractal phase diagram, whose landscape shows 
a qualitative agreement with experiments. 


The Fractional Quantum Hall Effect (FQHE) [T] is 
among the most fascinating quantum phenomena in¬ 
volving strongly correlated electrons. It attracts and 
fuels research in many directions since its discovery [5]. 
Lately, much interest has been directed to quantum Hall 
states as experimentally accessible prototypes of topo¬ 
logical states of matter, which have promising applica¬ 
tions to quantum computation I3H5]. 

The physics of the FQHE is well-understood phe¬ 
nomenologically thanks to the pioneering work by 
Laughlin and his celebrated ansatz for 1/m filling frac¬ 
tions [5] . The approach was generalized to more compli¬ 
cated fractions through the introduction of composite 
fermions [a [8] and a hierarchy of quasi-particles with 
fractional statistics [SHE], or by conformal invariance 
arguments [IMS]- A huge amount of results were ob¬ 
tained in the years, confirming the validity of the ap¬ 
proach based on model wavefunctions HdlHIS]- 

There is an ongoing effort toward the formulation of a 
systematic microscopic theory of the fractional quantum 
Hall effect. An intrinsic difficulty is the absence of an 
evident perturbative parameter, a common hindrance in 
strongly-correlated systems [5|. In 1983 Tao and Thou- 
less (TT) observed [20] that electrons in a strong mag¬ 
netic field could form a one-dimensional Wigner crys¬ 
tal m in the lattice of degenerate states in the lowest 
Landau level (LL), and suggested that this mechanism 
may explain the fractional quantization of the Hall resis¬ 
tivity. However, the resulting many-body ground state 
displays long-range spatial correlations, in conflict with 
Laughlin’s results. This route to a microscopic theory 
of the FQHE was abandoned (by Thouless himself [2^1. 
as the Laughlin ansatz offers several advantages, e.g. its 
high overlap with the exact low-density ground state, 
and the fact that it constrains very naturally the filling 
fractions to have odd denominators. The TT frame¬ 


work was recently reconsidered by Bergholtz and co¬ 
workers p5I125j . They found that TT states become 
the exact wavefunctions of the problem in the quasi one¬ 
dimensional (thin-torus) limit. 

Nowadays experiments in ultrahigh mobility 2D elec¬ 
tron systems are revealing a fractal scenario for the Hall 
resistivity as a function of the magnetic field: indeed 
more than fifty filling fractions are observed only in the 
lowest LL [26]. 

Here we study the thin-torus limit of the quantum 
Hall Hamiltonian in the lowest LL, and show that it 
realises a repulsive gas on the lattice of degenerate Lan¬ 
dau states, with the magnetic field acting as a chemical 
potential. The zero-temperature statistical mechanics 
of this class of models was studied extensively [271130] . 
It is characterized by an infinite series of second-order 
phase transitions, occurring at critical (non-universal) 
values of the chemical potential /i. The density of par¬ 
ticles p{pL) is the order parameter, and takes a different 
rational value in each phase, thus producing a devil’s 
staircase (a self-similar function with plateaux at ra¬ 
tional values also known as the Cantor function) when 
plotted against /i [55]. There is a revived interest in 
these models, for potential applications to quantum sim¬ 
ulators with ultracold Rydberg gases [3TH33] . 

Our mapping allows to (i) interpret the dependence 
of the transverse conductivity on the magnetic field as a 
fractal sequence of phase transitions, peculiar to ID re¬ 
pulsive lattice gases; (ii) establish the incompressibility 
of the ground-state hierarchy in the thin torus limit; (iii) 
provide a theoretical prediction of the relative widths of 
different Hall plateaux. 

We consider the standard two-dimensional gas of 
interacting electrons in a uniform positive background, 
providing charge neutrality. We make the assumptions 
that in strong magnetic fields the mixing between dif- 
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ferent Landau levels is suppressed, i.e. we work in the 
regime ^ tOc, where i = l/(ei?)^/^ is the mag¬ 
netic length and Wc = eBI'm is the cyclotron frequency 
(?i = c = 1) and spin degrees of freedom are frozen 
in the lowest spin level. We take the system to have 
area L^Ly and to be periodic in the y direction, so that 
the single-particle wave functions may be written in the 
form 


<Ps{x,y) 





( 1 ) 


with s = 1, 2,..., W = ■ The filling fraction v = 

Ne/Ns is less than one. 


In second quantisation, the Coulomb interaction be¬ 
tween the electrons in the lowest LL is 

Ns 

He— ^ ' ^SI —S3,S2—S3®Si®S2®®1+'52—S3®S3 ’ (^) 

31,S2,S3 = 1 

where aj, Os are fermionic creation and annihilation op¬ 
erators, and momentum conservation in the periodic di¬ 
rection is manifest. The Coulomb matrix element can 
be parametrized in a useful form by considering periodic 
boundary conditions in both directions (torus geometry) 
EOllMlES]- See also the Supplementary Material (SM). 




Si—S 3 ,S 2 —S 3 



The starting point of our analysis is the observation 
that this matrix element depends on a single variable 
in the thin-torus limit ijLx >> 1: the calculation (de¬ 
tailed in the SM) shows that the matrix element, when 
it is non zero, reduces to Vsi-S 3 ,s 2 -S 3 = e^/^W 4 i-s 3 
(with positive). By plugging this result into the 

Coulomb Hamiltonian we obtain 


He 


i 


^ Wsal^+,al^_,as^as, . 

Si,S2,S 


( 4 ) 


In the grand-canonical ensemble, the total Hamiltonian 
is the sum of the Coulomb term, the constant kinetic 
term and a term with chemical potential fl: 


S = 1 Si,S2,S 

( 5 ) 

where the definition /r(H) = (/i — Wc) highlights the 
dependence of the effective chemical potential on the 
magnetic field. Electrons in the lowest LL form a one 
dimensional lattice (that we call target space). Impor¬ 
tantly, they interact through a translational invariant 
interaction (in the target space). The Hamiltonian is 
diagonalized in the Fourier basis, where the creation op¬ 
erator for the mode k is c]. = e • We 

obtain the following diagonal Hamiltonian with periodic 


(^2 + ^ 


27rig^^(s2-S3) 

Ly 


4'7r2(3l-S3)2 

LI 


( 3 ) 


boundary conditions: 

-ffiLL =-At(H) ^ Tifc-I-y ^ w(\ki - k2\)'nk^nk2 , 
k—1 k\^k2 

( 6 ) 

with rik = c\ck and W{k) = e W{s) a repul¬ 
sive potential. The explicit form of IT(fc) is given in the 
SM; it decays as L^/{lk). 

This form of the Hamiltonian realises a mapping (in 
the thin torus limit L^/l 1) of the FQHE on a 
one-dimensional lattice gas with repulsive interactions, 
whose degrees of freedom are the Fourier modes of the 
target space. Notice that a generic quantum Hall Hamil¬ 
tonian on the torus is dual with respect to the unitary 
transformation defined by the Fourier modes, provided 
that Lx and Ly are exchanged (see the SM). In this re¬ 
spect our thin torus limit is equivalent to the one usually 
considered in the literature. 

As noted above, in these models the density as a func¬ 
tion of the chemical potential exhibits a devil’s staircase 
structure. Inspection of the Hamiltonian Q shows that 
the role of the density is played by the filling fraction 

whereas the chemical potential can be tuned by the 
magnetic field B. 

Schematically, the investigation of this class of mod¬ 
els follows two steps: (i) The ground state of the sys¬ 
tem is sought at fixed v = pfq {p and q coprime); this 
problem was solved by Hubbard [36]. (ii) The stability 
region A/i (under single particle/hole exchange) of each 
ground state is determined; this was done by Bak and 
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Bruinsma [55] , and by Burkov and Sinai [3D]. Both steps 
are subject to the technical condition that the potential 
be convex, which is fulfilled by the thin-torus potential 
W{k). We reproduce this two-step construction in the 
following. 

Intuitively, the ground state of a repulsive lattice gas 
at filling fraction v = p/q is, a, configuration where par¬ 
ticles are placed as far as possible from each other. The 
underlying lattice structure introduces the possibility of 
frustration, exhibited by deviations from the continuum 
equilibrium positions. The pattern of occupation num¬ 
bers can be obtained through the continued-fraction ex¬ 
pansion of v = p/q: 


no + 


’ • H- 

nx 

Each level in the expansion realises a better approxima¬ 
tion of I/; for rational v the number of levels A-l-1 is finite. 
At A = 0 (i.e. p =1), the ground state is a periodic crys¬ 
tal with inter-particle distance no = q, corresponding to 
Laughlin-type states. At A = 1 the inter-particle dis¬ 
tances can not be all equal, and a “defect” appears: the 
periodic ground state is formed by (ni — 1) Laughlin- 
type blocks of density 1/no and one block with density 
l/(no + 1); these correspond to Jain-type states (a con¬ 
cise representation is (no)"^“^(no -I- 1)). This construc¬ 
tion can be generalized iteratively to the level A (see 
Fig. [^for three examples, and the SM): the general rule 
uses the ground states at one level as building blocks to 
construct the ground states at the next level. The posi¬ 
tion of the j-th particle in the u = p/q ground state can 
be expressed compactly as lq/pj\ , where [-J denotes the 
integer part. (We notice en passant the connection with 
the sequences of characters known as Sturmian words.) 

Due to the periodic boundary conditions, the ground 
state at filling factor v = p/q has a g-fold degeneracy, 
corresponding to the possible translations in the target 
space. This plays an important role when quantum ef¬ 
fects are taken into account (see below). Summing up 
the foregoing observations, a compact form of our wave 
functions is the following: 

lpN„/q\ 

\n=p/q)r= n r = Q,---,q-l. 

i=i 

( 8 ) 

Once the ground states at general v have been de¬ 
termined, their stability under single particle/hole ex¬ 
change can be established. The stability interval in the 


V Cy configuration 
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FIG. 1. Hubbard ground states for different filling frac¬ 
tions V and their explicit periodic structure. The first two 
from the top belong respectively to Laughlin {v — 1/3) and 
Jain {v = 3/7) series. Each periodic configuration may be 
expressed in a compact way through the sequence Cv of its 
interparticle distances (the general algorithm to construct 
Hubbard ground states is extensively reviewed in the SM). 

effective chemical potential is given by EDI 

oo 

Ap{p/q) = 2q ^ k ^IT {qk + \) +W {qk — 1) — 2W (qk)'^ . 

k=l 

(9) 

As A/r.(j/) > 0 for all rational filling fractions, this con¬ 
struction yields a phase diagram where each rational v 
appears as the stable density for a finite interval of p 
(hence of B), thus realizing a devil’s staircase. As a 
consequence of our mapping, the stability equation (§ 
constitutes a proof of the incompressibility of the hier¬ 
archical ground states obtained in the thin torus limit. 

It is worth remarking that the precise form of the po¬ 
tential does not affect qualitatively this result, as far as 
the convexity condition is fulfilled. 

Our results support a new interpretation of the FQHE 
landscape (at least in the thin torus limit) as the 
zero-temperature phase diagram of a fermionic one¬ 
dimensional lattice gas model with repulsive interac¬ 
tions. The results reported above allow to plot a snap¬ 
shot of the relation between magnetic field and inverse 
filling fraction. To this end, we assume that even- 
denominator ground states, which are not seen in the ex¬ 
periments, are gapless. A possible argument, related to 
the magnetic translation group symmetry, has been pro¬ 
posed by Seidel m (see the SM). With this assumption, 
we set Ap{p/q) = 0 for even g, and use the stability for¬ 
mula ^ otherwise. The potential W has a non-trivial 
dependence on the magnetic length As noted above, 
it decays algebraically as l/{^k). To obtain a large dis¬ 
tance ^-independent behaviour, the chemical potential 
is rescaled as p ^ pi"^ , which is equivalent to a rescaling 
of the entire Hamiltonian, H —)■ Operatively, we 

set a cutoff gmax on the possible denominators, we list 
(in increasing order) all filling fractions p/q such that 
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magnetic field [arbitrary units] 


FIG. 2. Inverse filling fraction l/v plotted against the mag¬ 
netic field B (in arbitrary units). The most visible plateaux 
are highlighted with their corresponding occupational peri¬ 
odic pattern in the reciprocal target space. This snapshot 
shows a qualitative agreement with the experimental mea¬ 
sures of Hall resistivity, both for the relative widths of the 
plateaux and for the quasi linear trend of the landscape as 
a function of B. In the inset, a portion of the staircase is 
magnified and some experimentally-observed plateaux 
are marked. 


q is odd, q < qmax, and p = 1,..., g, and we compute 
A/i for each one of them. Doing this by increasing order 
allows to obtain iteratively the two stability boundaries 
P- and of each plateau; the corresponding values of 
the magnetic field and B+ are calculated from the 
relation p, = —p/{eB) — 1/m. The resulting landscape, 
presented in Fig. is qualitatively in accord with the 
well-known behavior obtained in experiments. 

The roles of the numerators and the denominators 
in the filling fractions have competing effects on the 
plateau widths. Equation implies that the width 
of a plateau (in the chemical potential p) only depends 
on the denominator. Filling fractions with the same de¬ 
nominator will have the same lS.p. In particular, it can 
be easily shown by use of Eq. 0 that the plateaux get 
narrower as the denominator q is increased. However, 
the non linear dependence of ^ on i? breaks this sym¬ 
metry, by enhancing the stability of plateaux at larger 
magnetic fields. As a consequence, filling fractions with 
the same denominator have larger stability intervals (in 
B) for smaller numerators p. The most evident exam¬ 
ple of this general mechanism can be recognized in the 
fact that the plateau at jz = 1/3 is larger than that at 
^ = 2/3, as is experimentally observed. 

Notice that, in statistical mechanics, systems with 
slowly decaying potentials are pathological: their free 
energy is not extensive as a function of the particle 


number. In our framework, this has the effect to push 
the staircase toward infinity as the cutoff gmax is in¬ 
creased. This issue may be overcome by regularizing the 
Coulomb potential. Our thin torus analysis is largely 
independent of the precise form of the potential. 

We remark that the continued-fraction expansion that 
we employ to construct the ground states naturally pro¬ 
vides a definition of “complexity” of a given filling frac¬ 
tion, via its level A. This construction has a natural 
interpretation in terms of quasi-particles [38] . that we 
do not further pursue here. 

The main result of this work is the mapping between 
the Hall Hamiltonian in the thin-torus limit and a long- 
range repulsive lattice gas model in one dimension. This 
results allows us to interpret the FQH ground states 
as Hubbard states, and to prove their incompressibil¬ 
ity, as a direct consequence of Eq. §. The lattice gas 
also brings to a scenario where the Hall resistivity as a 
function of the magnetic field is a devil’s staircase. By 
assuming that even-denominator ground states are gap¬ 
less, qualitative accordance with the experimental land¬ 
scape is obtained. This suggests that it may be fruit¬ 
ful to investigate the nature of the correlated ground 
states at more exotic fillings in the lowest LL. This 
is in principle possible by generalizing the composite- 
fermion picture (recently used to propose new incom¬ 
pressible ground states at v = 4/II and v = 5/13 [3^). 
or by exploiting the recent results with Jack polynomi¬ 
als [T^liOlliT] . 

We are grateful to Bruno Bassetti, Sergio Caracciolo, 
Mario Raciti, Marco Cosentino Lagomarsino, Andrea 
Sportiello, and Alessio Celi for useful discussions and 
advice. 
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I. DERIVATION OF THE YOSHIOKA AND TAO-THOULESS FORMULAE FOR THE LLL COULOMB 

MATRIX ELEMENT 

The earliest numerical calculations of the ground state for FQHE including Coulomb interactions for various filling 
fractions were done by Yoshioka et al. [1, 2]. Here we provide the details of their evaluation of eq.(2), in the lowest 
Landau level. The geometry is that of a periodic array of rectangles with sides and Ly and area L^Ly = 2Tr£‘^Ns, 
where Ns is a natural number; uj is the rectangle [0, x [0, Ly]. From the Ns degenerate LLL eigenstates with 
centers in w 


i’sir) 




1 . 

( X 

27rf ^ 

27r 

~2 ^ 


■ ~r^ 

Ay y 

-i — sy 

l^y 


0 < s < Ns - 1 


( 1 ) 


one constructs a basis of orthonormal quasi-periodic eigenstates: 0s(r) = (r). 

An electron in oj interacts with electrons in oj as well as with their copies. The Coulomb interaction depends on 
r = ri — r 2 and is a periodic function of the lattice; 


;(r) = 


E 


+ nixLx)'^ + {y + myLyY 


It has Fourier expansion v(r) = 77 ^ where Qx = and Qy = ^Uy and u(q) = f^drv(r)e = 

n 2 ^ ^ 

. In the Fourier representation the integrals for Coulomb matrix elements factorise: 


(S1,S2|U|S3,S4) 



dridr2 6»si(ri) 9s^{r2)v{ri 


r2)0ssiri)0sA^2) 


-j^ E ^'(q)-^i.3(q)^2.4(-q) 


The integrals Is^s' (q) are independent of the potential, and are now evaluated: 


Is, s' (q) = / dr0s {r)0s' (r) exp(zq • r) = e 

J (jj 

where 5' means equality modulo Ns. 




—s'+ny,0 
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Proof. The integral in y is straightforward: 


^ 6s- 
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/■ “= dx -i [2 
s' + {m—m^)Ns+ny^0 I n t —^ 
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sY-^s'Y +iq. 
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^ (s + S ) 2 *?a: 


The double sum involves m+m' = y and m—m' = v, and f{m+m', m—m') = ^ /(2/i, 2:/)+/(2/i+l, 2i^+l). 

Therefore: 


7^,s/(q) =e 




'^T[' 

/7,l/ 


/•ix 
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Os-s'+2i/Afs+ny,0 / . L y J 


+ ds-s' + (2i/+l)Ar,+ny,0 ^ ^ 


The sum on produces a Gaussian integral on M. The two integrals have the same value. The final result is 
obtained. □ 


The Coulomb matrix element is: 


(Si, S2|U|S3, S4) = -j—j- X] 


^xLy |q| 


„-5lq|^^^+*9x^(si+S3-S2-S4) r/ T! 

^si —sa+riy ,0^S2 —S4 —ny ,0 


The two constraints imply momentum conservation: si + S 2 = S 3 + S 4 modulo fds- Eq. (2.9) in Yoshioka’s paper [2] 
is obtained: 


(si, S2|X|S3, S4) = V 


LxLy ^ |q| 


S3 — Sl,qyLyj2TT 


( 2 ) 


It is an exact formula. The Tao-Thouless formula is now obtained. First use the constraint 6' to sum on 

-^^ffiss-si+mNsp 


(S1,S2|X|S3,S4) = ^Sl+S2.S3 + S4 X T 6 


- V9x+®9x^£f-(«3-S2) 


LxLy 
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+ ^{S'i- Si+ mNs^ 


Next, approximate the sum on by an integral ~ neglect terms m ^ 0 because of the exp factor. 

Eq.(3) in Tao and Thouless, [3] is obtained: 


g2 _2lL^(g2_g^)2 POO g“T9 +i9T^(33-S2 

(si,S2|u|s3,S4) = —(531+82.33+346 / dq 
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+ ^(S3 - Sl)2 


( 3 ) 


With q£ = t the matrix element is (conservation of momentum is not specified): 


2e2 -^( 33 - 31 ) 

V = —e y 
J-Jt, 


dt 


e 2 


: COS t 




^(63 -Sl) 


Lx S 3 — S 2 


n 2 \ e Ns 


II. THIN TORUS LIMIT 


In the limit Ly, <C 7 (thin torus limit) the cosine function equals one. The integral is a Bessel function 

E(si 3) = f^(s3 - Si)" 

Ly V L^ 
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For the largest values of S 13 , the argument of Kq is of order <C 1. Therefore, only the small argument behaviour 

matters: 


-y- 

J-jq, 


( or 2 ^13 ) + C' + vanishing terms 

V 2T 


(4) 


where C is Euler’s constant. The Coulomb term of the second-quantised Hamiltonian in the thin torus limit is 
diagonalised by a change of basis: 


iVs 
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V=4 y] F(si3)a^^a^^Osi+s2-s3ds3 = I ^ Vk-k'nkfik'- ^VqN 
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Si,<52 ,53 = 1 
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where fik = c\ck, 



( 6 ) 


and Vk = / 2 ^(®) exp(i|pfcs). For large Ns, the Fourier series is approximated by an integral (s = xNg, 

1 = Ngdx): 
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Jo dx log X cosiirkx) = 4 dt^, k^O 

The Sine Integral function slightly deviates from it asymptotic value tt/2 only for fc < 1. We then approximate: 
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III. SYMMETRIES OF THE EXACT LLL HAMILTONIAN 

J_ ^ r°° du 


In the exact formula for the matrix element by Yoshioka, (2) insert = f -4^e 

J y ■/ |q|£ J-00 


T 1234 = 
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Sl—S 3 , —"y S2—S4 ,"h 


Since the delta-contraints and the exponentials involve nx,y modulo Ns, let Ux = RNs + ^, Uy = SNs + v, R, S G Z, 
fi,!/ = 0 ... Ns — 1. The point {nx,ny) = (0, 0) i.e. i 2 = /r = i? = S' = 0is forbidden. 


Ns-l 




du 

— 00 's/ ^TT 


00 00 


/ ft n n 


^-nN,^(R+^)Rl+u^)-^N,^iS+Tyfil+u^)j2nt.^c> 


Si 


Si—S 3 , —S 4 ,i^ 


R— — (yc) S— — (X) 


In the sum, let’s single out fj, = 12 = 0 (then RS 7 ^ 0) from the values where fi and v are not simultaneously zero. The 
sums can be expressed as Jacobi Theta functions: 


V1234 = 


du 1 


v^i + uH 

Ns-l 

E 

fi.i'—O 


^0 0 


iLx 1 


LyNs 1 + 


d3(0' ^ 


-1 


^Si—S3,0^S2—S4,0 



iLx 1 

^3 

f HU 

iLy 1 \ 

\Ns 

LyNsl + u'^J 

\Ns 

LxNsl+u-^J 


LxNg \ 


'^Si—S3, —IV^S2—S4,l/ 


Ns-l 


— 53 ,0^S2 —S 4,0 ' / ^ ^ \t^'> ^ S4,L 

fi,y—0 
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where: 


Vb = 


1 


dt 




2N^V^ Jo V^(l-t) L 

dt 


^3 0 


iLxt 


LyNg 


^3 ( 01 


LrNg 


- 1 


1 


2N^V^Jo y^t{l - t) 


/ TT^ 

iLxt 

^3 

(-nv 

iLyt \ 

\Ng 

LyNg) 

KNs 

LxNg) 


with the periodicities V{y, + Ng,v) = V{Ng — ^,u) = V{fi, v), and the same for v. 
In second quantisation the Coulomb operator is: 


Ns 


V = i 
2 


^ (si,S2|u|s3,S4)a^^a^2as^as3 

SiS2<S3'54 = 1 

2 2 ^ s — ^ ^ s 

fj.,iy=0 Sl234 = l 


( 8 ) 

(9) 


V=^Co(iV^ 


2 Ns-1 Ns 

S2,S3 = 1 


( 10 ) 


An alternative expression arises by representing the delta functions as Fourier sums: 

2 Ns-l Ns 


v=5jr„(iv"-jv) + - v{„.s) Y. o 

— 0 Sl234 = l 


12-77fl - 


1 

ivF 


Ns 

E 

p,k—l 


J^k{si-S3 + iy)+i^p(s2-S4-iy) i t 


It is natural to introduce the new canonical basis of operators (6), fc = 0,..., — 1, 


V = 


2 Ns-l Ns 

p,k—l 


( 11 ) 


The expressions (10) and (11) are exact representations of the Coulomb operator in two basis related by the discrete 
Fourier transform (6). They differ in the basis operators and the exchange of /r with v, i.e. between and Ly. 

As the two canonical basis are related by a unitary transformation, 


uV,u = cI 


it is 


UlH(r)U = Hf- 

where H(r) = K + V(r), K = IfiwcN, N = 4afc. 

Two unitary groups of translations are introduced: 


T = exp 

They have actions (modulo Ng): 


/carafe 


fc=0 


S = exp 


Ns 


1E^4c. 




rp t rpt _ J 


= 4+1 


T4Tl = e*^'=4. SlclS = e 




TU = US 


( 12 ) 


(13) 

(14) 
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and fulfil the following algebra, 



Because of momentum conservation, it is: 



for any aspect ratio r. Therefore the eigenspaces of H(r) are left invariant by both unitary operators. 


(15) 

(16) 


Proposition 

Let H be a Hamiltonian that commutes with T, S, N, and S be the eigenspace of ground states of H with N particles, 
with projector P. v = p/q and £ has exactly dimension q, then: 

1) PT«P = PAt, PS«P = PXs, |Ar| = |As| = 1; 

2) If PUP ^ 0 then | At = As ■ 

3) If [H, U^] = 0 then At = ±1 and As = ±1 

Proof. The operators H,N and T can be diagonalized simultaneously, with eigenvectors Because of (15), 

the q vectors t), r = 0,..., g — 1, are orthogonal eigenvectors of T with eigenvalues E, v and t exp(—z27rri^). 

Therefore they span £. 

1) It follows that the operators T'^ and act on £ as multiples of the identity, with eigenvalues At and As- 

2) From the relation U^T^U = S« one gets: PT«UP = PUS«P. Therefore: 

(At-A s) PUP = 0 

If PUP ^ 0 then the eigenvalues are equal. 

3) If [H, U^] = 0, the set £ is left invariant by U^. Since U^T'^U^ = it is At = ±1. The same is true for As- □ 


IV. SEIDEL ARGUMENT FOR ODD DENOMINATORS 

The proposition applies to the Hamiltonian H(r) for any r. It becomes more stringent in the thin-torus limits 
r —)■ 0 and r —)■ oo, as discussed by Seidel [4] for discarding filling fractions with even denominators. 

In the thin torus limit, the eigenvalue problem is mapped to a lattice gas eigenvalue problem that was analytically 
solved by Hubbard. The eigenspaces of the unitarily equivalent Hamiltonians H(0) and H(oo) consist of occupation 
number vectors in the basis a\ and cj, respectively. 

At filling fraction v = p/q, the ground states |ni, n2,..., u-tv^) are characterized by periodic sequences rir+q = rir, 
with ni + • • • + Uq = p, Ng = Mq and N = pM. For the two Hamiltonians the numbers are the same, but refer to 
occupations of basis vectors related by a Fourier transform. 

On a ground state of H(0), \n) = |0), the operator T is diagonal, while S acts by shifting 

by one all numbers, giving a new ground state. As a result of Hubbard’s prescription, the eigenspace £(0) (and i?(oo)) 
has dimension q. The operators PT'^P and PS'^P are multiple of the identity: 

0 _ “-1 9 

PT«P|n)=exp iq^^^{jq + k)nu |n) = (-1)‘^^^^"^V) (17) 

® j=0 fc=l 

PS«P|n) = al”\ ... 4”’ ... 4”^ ... 4”’ |0) = (-1)P'(“-1) \n) (fermions) (18) 

By the previous proposition, if PUP 0, it is (—= (—1 )p This is verified if p is even, since q is 

forced to be odd (they are coprime). If p is odd it must be (—= (—Since M is a scale factor, the 
statement independent of it is that q must be ODD (for fermions). 






6 


V. EXPLICIT CONSTRUCTION OF THE HUBBARD GROUND STATES 

In a remarkable paper Hubbard [5] studied orderings of electrons in quasi-one-dimensional conductors. When the 
electrostatic energy is dominant, a ID lattice-gas model results, and Hubbard gave a general method to determine 
the exact ground state. 

Consider a chain of N sites. The particle configuration is a vector n of occupation numbers = 0,1 (t = 1,..., N). 
Particles interact via two-body forces that depend on their distance, and not with themselves. Hubbard’s Hamiltonian 
is: 

= j\)n,nj (19) 

3¥^i 

We search for the configuration (GS) that minimizes E(n), at fixed number of particles m = or density 

p := m/N. Since the potential is repulsive, intuition suggests to allocate particles on the lattice as far as possible from 
each other. However, the prescription of the density does not allow to do this evenly. Hubbard analitically solved the 
problem for the infinite chain and for potentials satisfying the two conditions: 

U(r)—as r —>■ oo (20) 

U(r-I-1)-I-U(r — 1) > 2U(r) for all r>l (21) 

The convexity condition is satisfied by the Coulomb potential and by V{\i — j\) = |i— j|““, a > 0. Hubbard’s solutions 
are called generalized Wigner lattices or most uniform configurations. They are independent of any further detail of 
the interaction potential. 

Let’s extend the A-site chain into an infinite one with particle density p, and search for ground states that are 
fV-periodic, with pN particles per period. We thus impose the periodic boundary conditions (pbc) 

^ — f j ■ ■ ■ 5 (^^) 


and the bound 

N+j 

'^n,=m Vj 

i^j 

(23) 

The A-site chain is a loop, and it interacts with the rest of the lattice, i.e. with an infinite number of self-copies. A 
complete rotation of the loop corresponds to a shift of one period of the infinite chain. 

Let’s label the particles round the loop as ^ = 1, ... ,m and define r^'^ as the position of particle v, and ri^'' as the 
interval between particles ly and ^ -k 1. Then ri^^ =: r^J^h — ri^K Because of pbc, we can obtain the position of a 

generic particle in the infinite chain, i.e. -k N. It follows that = ri^K 

The distance between particle ly and v + p is 

„(u) — ^(0) 

' V ' ' 

T — U T—0 

(24) 

From this we obtain the useful properties 



ri^+m = &'> for all M 

(25) 

(26) 


The energy only depends on relative positions, {rl^^}. Therefore the vector =: {r^^} defines a system’s 
configuration and satisfies defines a configuration up to transformations like ri}^ —)• and a 

configuration satisfies implies 

How can we express energy (19) in term of The interaction energy between particle v and particle iz + p, is 

V{rlf''^). Therefore, the interaction energy between particle n and the rest is 


Viri^^) + U(r^''+™)) -k U(r(,^+2™)) -k ... 

fi fj,—! k fi—1 
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{ui} configuration 


1111100000000 • 

• 

• 

• 

• 

O 

O 

O 

O 

O 

O 

O 

o 

{1,1,1,1,9} 

1001001001010 • 

o 

0 

• 

O 

O 

• 

O 

O 

• 

O 

• 

o 

(3,3,3,2,21 

0000110100110 o 

o 

0 

o 

• 

• 

O 

• 

O 

O 

• 

• 

o 

(1,2,3,1,6} 

1010010100100 • 

o 

• 

0 

O 

• 

O 

• 

O 

O 

• 

o 

o 

12,3,2,3,3} 


TABLE I: Examples of microscopic configurations for a loop of A = 13 sites, p = 5/13 


that we can simplify using property (26) 


= +kN) (27) 

k 


The energy of the system 

1/ fl 

^ m m 

= 2 E E E + kN) + l/(ryy + kN) + ... 

v—1 k fi—1 

can be simplified using property (25) 


E{r<^^'>) = E E +kN)=AY,Yl (28) 

k i^—l k 

where A it is a factor that depends on lenght of the chain, =: {T^} and 

m 

[/fc(r('')) =: y(r(,^) + kN) (29) 

V=1 

We now have the basic ingredients to solve the problem. The first task is to look for an algorithm to minimize (28) 
with the condition 


^rW=iV (30) 

This problem is similar to m problems in which we separately minimize the inner sum (29) with the condition 

m 

Y&^=^^N (31) 

1^=1 

If such a solution exists, it also minimizes (28). To obtain the solution for inner problems we need the following 
theorem: 

Theorem V.l. If {vi} = ri,..., is a set of m integers such that 

m 

■= R = mr + a, 0<a<m (32) 

i=l 

and if V : N 1 -^ h is an integer function such that it is strictly convex, then 

m 

(m — a)V (r) + aV (r + 1) < 

2 = 1 


(33) 




Proof. First we suppose that {rj} is a set such that — rj\ < 1 for all pairs i,j; such a set will be called minimal. 
We define f =: min(ri). For a minimal set we have f < < f + 1 for all v. can assume only values f or f + 1. Let 

h be the number of that assume the value r. Then m — h is the number of that assume value f + 1. Calculating 

m 

Ti = nf + {m — n){f + 1) = mr + m — n (34) 

Z = 1 


by the uniqueness of the modular decomposition, we obtain r = r and m — n = a. So it follows at once that, for 
minimal sets, the equality in (33) is satisfied. 

Now let C be a non minimal set {ri}; then, for some s 7 ^ t one has rg > r* + 1. We construct the set C by taking 
C and moving only Vg and rt'. r' =: for alH 7 ^ s,t, r(, =: — 1, =: rt + 1. We see that C satisfy the condition 

(32) and, by using the lemma (V.2), we obtain 

m m 

V{r[)pV{r',)<V{rg) + V{rt) ^ C(r') < ^ C(r.) (35) 

i=l i=l 

If C isn’t minimal we repeat the procedure to obtain a C" with 

until one arrives at a minimal set (7°. Using inequalities between ,... ,C one proves the thesis □ 


Lemma V.2. Let V : N 1 —)■ Z on integer function such that it is strictly convex i.e. U(r + 1) + U(r — 1) > 
2 U(r) for all r > 1 then 

U(s + 1) + U(t — 1) < U(s) + U(t) for all s such that s <t (36) 

%endalign* 

If it is a minimal set, it minimizes the energy (29). We need to fabricate such set. Starting from (31), we 

calculate integers and by decomposing the fraction /i/p into the sum of a non-zero integer and a proper 
fraction: 


E 


= mr^^^ + 


p m 


(37) 


We obtain r^x-) = [/r/pj ([xj means interger part of x). If ri^^ e S^{p) := ^r^t^\r^x-) _(- ij = [[/i/pj, [p/pj + 1] then 

{ri^^} it is a minimal set. There are m — of ri^'^ that take value [p/pj and of ri^^ that take value [p/pj + 1. 
Therefore a set of solutions that minimize (29) exists for each p = 1,... ,to, and for all fixed p, the energy (28) is 
minimised when r[P € S^{p) for each p = 1,..., m. 

Example V.l. to = 5 and N = 13 (p = 5/13/. 

Calculating r^^^ and we obtain minimal configuration for each p. Results are in the table. 


h 

j.(m) 

TO — 

j.(v) + 1 



1 

2 

2 

3 

3 

{ 2 , 2 , 3 , 3 , 3 } 

2 

5 

4 

6 

1 

{ 5 , 5 , 5 , 5 , 6 } 

3 

7 

1 

8 

4 

{ 7 , 8 , 8 , 8 , 8 } 

4 

10 

3 

11 

2 

{ 10 , 10 , 10 , 11 , 11 } 

5 

13 

5 

14 

0 

{ 13 , 13 , 13 , 13 , 13 } 


TABLE II: Minimal sets for U, m — 5 and N = 13 

Each set of r^^ minimizes the energy (29) indipendently from the internal order of numbers. For example both 
{2,2,3,3, 3} and {2,3, 2, 3, 3} minimize U{r^^^). But from {r[}^} = (2,2, 3, 3,3} it descends that the configu¬ 

ration {4,5,6, 6 , 5} is not a minimal set for I7(r*^^^). Therefore, to minimize energy (28) we need to combine results 
in table and choose only one {rl^^} configuration, i.e. to .specify the ri^^ ordering, up to rotations. 

In the example we have only two ineguivalent configurations. They generate the chain of configurations 

{2,2,3,3,3} ^ {4,5, 6 , 6 ,5} ^ {7, 8 ,9, 8 ,7} ^ {10,11,11,10,10} ^ {13,13,13,13,13} 

{2,3,2,3,3} ^ {5,5,5, 6 ,5} ^ {7, 8 , 8 , 8 , 8 } ^ {10,11,10,11,10} ^ {13,13,13,13,13} 

The configuration that minimises the total energy is {rl^^} = {2,3, 2, 3, 3}. It is the most uniform configuration for 
p = 5/13. We can write this result in a compact notation: {2, 3, 2,3,3} =: (10100)^(100) = (23)^(3) 



9 


We proved that, for each /r there is a set of configurations that minimize and we know how to construct it. 

We don’t have yet a method to choose the one that minimises the total energy (28), but we can fabricate it. We start 
using /i = 1, i.e. we minimize the nearest-neightbor energy We obtain the set of minimal configurations by 

applying the algorithm of the theorem. We calculate and from the decomposition 


1 

P 


= r(i) 


m 


(38) 


and then we have the set of configurations that minimises i.e. such that there are m — times 

and times + 1. 

For example, if = 0, i.e. p = (the simplest possibility), we have only one configuration that satisfies the 

problem: particles are distributed uniformly, one every sites. Indeed ..., This 

microstate solves the general problem, as we can verify by comparing the chain of configurations that the solution 
generates with sets of configurations that minimize U{r^^^): 


... ,r«} ^ {r^^) = ... ,r(2)} 


Let’s introduce a new notation to indicate the set of configurations that minimises U{r^^''). An element can be 
j-(i) or + 1. To identify only one system’s configuration we have to know the pattern of in the vector 
up to rotations. However, knowing that means that after particle v we must have — 1 holes. In terms 

of occupation numbers, when we find we fix a part of the configuration, that we indicate with round brackets. 
We have only two types of partial configurations: (10^'’^ *~^^) and (10^’’* *^). If we assume that a particle is always 
at the beginning of the string, we can compact the notation for partial configurations using (10^'’* *“^^) =: (r*-^^) 
or (lO*-’"'' ^) =: (r*-^^ + 1). A choice of configuration means fixing the order of partial configurations. A minimal 
configuration for U{r^^'>) is part of 


Ml 


{ni\ 


(r^^)) X [m — 1 

(7^(1) + 1) X J 


(39) 


To understand the notation let us consider examples. 

Example V.2. m = 5 and N = 13 {p = 5/13). We obtain = 2, so the partial configurations are (10) and 
(100). Minimal configurations for {/(r^^)) are {2, 3, 2,3,3} = (10)(100)(10)(100)(100) = (23)2(3) and {2,2,3, 3, 3} = 
( 10 )( 10 )( 100 )( 100 )( 100 ) = ( 2 ) 2 ( 3 ) 3 . 

Example V.3. p = 1/n {m = N/n). We obtain = n, so the partial configurations are (10"“^) and (10"). 
Minimal configuration for U(r^^'') is {n,... ,n} = (10"“^)’" = (n)™. 


Let us observe that the smallest possible partial configurations are (1) and (0) (that we call 0-particles). Then, a 
generic configuration is part of: 


- {<"•) ^ (?) f M } («' 

Knowing the minimal sets for U{r^^'^) (i.e. p density decomposition) permits a first ordering of 0-particles in bigger 
partial configurations Mi C Mg: 


{Ui} 


(0) X [A^ — to] 
(1) X [to] 



(r^^^) X [to — 1 

(7.(1) + 1) X [a^^^] J 


(41) 


How can we select the configuration that minimises (28) among those (39) that minimise [/(r^^^)? If = 0 we 
know the solution: we have only one configuration in the set (39) that is a GS. If 0 the second term /m 

in the continued fraction decomposition represents the density of (r*-^^ -I- l)-site intervals which can be regarded as 
particles, within the (r^^))-site intervals, which can be regarded as holes [6]. In this way we define two types of 
1-particles. Therefore we want to find Hubbard GS for a density p = /to. We know which is the first step to solve 
this problem: we calculate and b^^l from the decomposition 


TO 



-f 


5(1) 


( 1 ) 


TO 


(42) 
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If is the vector of distance between 1-particles, i.e. it indicates how many sites intervals there are between 

_|_ 1 ) sites intervals, we have m — of that are and of that are -I- 1. So we obtain two types 
of 2 -particles 




and a subset of Mi that minimises also 

Ma := 


r 1 -I-x [m — 1 

' ((rW + l)(r(i))(«‘^’))x [6W] j 


(43) 


(44) 


If = 0 we have only one conhguration that is a solution. If b^^'> ^ 0 the second term jm in the jm density 
decomposition represents the density of (1) 2-particles within the (0) 2-particles. We can repeat the algorithm used 
before for 1-particles. And so on, as long as we find a null density fraction for /c-particles. Indeed in such case we 
have only one configuration in M^, that is a GS. 

Example V.4. m = 5 and N = 13 {p = 5/13). We start from the set of 0-particles 

M - //„ r . (0) X [8] 1 

Mo .- |{n,| . ^ jgj I 

From 5/13 decomposition we obtain the set of 1-particles 


13 „ 3 

— = 2 -I- - 
5 5 


Mi:=Una : 


(2) X [2] 

(3) X [3] 


^ = 1 + ? 
3 3 


From 3/5 decomposition we obtain the set of 2-particles 

M 2 ■■= |{ni} 

From 2/3 decomposition we obtain the set of 3-particles 

I = 1 + i ^ M. := {{».) 

Finally from 1/2 decomposition we obtain the set of A-particles 

0 


(3) X [1] 

(32) X [2] 


(32) X [1] 
(323) X [1] 


2 = 2 +- 


M 4 := {{nj : (32332) x [1]} 
Example V.5. m = 11 and iV = 47 {p = 11/47). We start from the set of 0-particles 


Mq ■■= < {ui} : 


(0) X [36] 
(1) X [11] 


From 11/47 decomposition we obtain the set of 1-particles 


^=4+^-^M4:=UnJ 


. (4) X [8] 

■ (5) X [3] 


11 o 2 

y “^^3 


From 3/11 decomposition we obtain the set of 2-particles 

M 2 ■■= { {nj 

From 2/3 decomposition we obtain the set of 3-particles 

M 3 := / {m} : 


3 , 1 

2 “^+2 


Finally from 1/2 decomposition we obtain the set of A-particles 

2 = 2 + Y —s> 


(544) X [1] 
(5444) X [2] 


(5444) X [1] 
(5444544) x [1] 


M 4 := {{n^} : (54445445444) x [1]} 
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Remarks: 

1) The algorithm has an end. As we are representing the rational number p as a continued fraction, it is known that 
every finite continued fraction is a rational number, and a rational number can be represented by a finite continued 
fraction in precisely two different ways. 

2) If p = {mp)/{Np) where m and N are coprime, and p > 1, the GS that is obtained by the algorithm is the same 
one as p = p/N. Then it is sufficient to obtain GS for densities in which numerator and denominator are coprime. 

3) Particle-hole symmetry makes it sufficient to obtain the GS for p < ^ [6]. 


Let us summarize Hubbard’s algorithm for a filling fraction p. First p is represented as a continued fraction, 


P = 


1 


Ml + 


1 

ll2 + ... 


[0;ui, ...,Uk\ 


(45) 


Then we define n-particle Yn and n-hole A„ recursively 

(46) 

(47) 

The initial conditions are 0-particle and 0-hole, Aq := (0) and Yq := (1). Hubbard’s ground state is X/,. 


A„ := (y„_i)(A„_i)“"-i 
:= (y„_i)(A„_i)“" 


p Xk configuration 


1 (3) 

• 

0 

o 

• 

o 

o 

• 

o 

o 

• 

o 

0 

0 ... 

1 (322) 

• 

0 

o 

• 

o 

• 

o 

• 

o 

o 

• 

o 

• - - • 

X (32)2(3) 

• 

o 

o 

• 

o 

• 

o 

o 

• 

o 

• 

0 

0 • • • 


TABLE III: Hubbard ground states 
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